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ABSTRACT 

We present the analysis of Herschel SPIRE far-infrared (FIR) observations of the z = 2.515 lensed galaxy 
SMM J163554. 2+661225. Combining new 250, 350, and 500 /im observations with existing data, we make 
an improved fit to the FIR spectral energy distribution (SED) of this galaxy. We find a total infrared (IR) 
luminosity of L(8-1000 jim) = 6.9 ± 0.6 x 10" L©; a factor of 3 more precise over previous L/r estimates for 
this galaxy, and one of the most accurate measurements for any galaxy at these redshifts. This FIR luminosity 
implies an unlensed star formation rate (SFR) for this galaxy of 119 ± 10 M© yr"', which is a factor of 1.9 ± 
0.35 lower than the SFR derived from the nebular Paa emission line (a 2.5(j discrepancy). Both SFR indicators 
assume identical Salpeter initial mass functions (IMF) with slope F = 2.35 over a mass range of 0.1 - 100 M©, 
thus this discrepancy suggests that more ionizing photons may be necessary to account for the higher Paa- 
derived SFR. We examine a number of scenarios and find that the observations can be explained with a varying 
star formation history (SFH) due to an increasing star formation rate (SFR), paired with a slight flattening of 
the IMF. If the SFR is constant in time, then larger changes need to be made to the IMF by either increasing 
the upper-mass cutoff to ^200 M©, or a flattening of the IMF slope to 1.9 ± 0.15, or a combination of the 
two. These scenarios result in up to double the number of stars with masses above 20 M©, which produce the 
requisite increase in ionizing photons over a Salpeter IMF with a constant SFH. 

Subject headings: infrared: galaxies — galaxies: high-redshift — galaxies: starburst — galaxies: individual 
(SMM J163554.2+661225) 



1. INTRODUCTION 

The recent launch of the Herschel Space Observatory (Pil- 
bratt et al. 2010) has opened the high-redshift universe to de- 
tailed far-infrared (FIR) investigations. Herschel observations 
of galaxies are specifically useful as they constrain the peak of 
the FIR spectral energy distributions (SEDs) for z^2 galax- 
ies, allowing robust measurements of the total infrared (IR) 
luminosity, L/^ = L(8-1000 /im), and dust temperature (e.g., 
Amblard et al. 2010; Elbaz et al. 2010; Magdis et al. 2010). 
From the total IR luminosity one can infer the star formation 
rate (SFR) and compare to SFRs based on other indicators. 
Both L/R and the dust temperature allow for detailed physical 
investigations of star-forming galaxies not possible without 
FIR observations. 

One important factor for understanding star-forming galax- 
ies at high-redshift is the form of the initial mass function 
(IMF), specifically at the high-mass end. There is currently 
an ongoing debate as to whether the IMF is universal. In 
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nearby galaxies the IMF is well fit by a Salpeter (1955) IMF, 
where dN/dM oc M"^, with slope F = 2.35 (Elmegreen 2006; 
and references therein), or even steeper (F > 2.35) in a few 
examples, such as in regions of the Large Magellanic Cloud 
(Gouliermis et al. 2005) and low surface brightness galaxies 
(Lee et al. 2004). On the other hand, the summed IMF in clus- 
ters of galaxies appears to be Salpeter or slightly flatter (F < 
2.35) than Salpeter (e.g., Renzini et al. 1993; Lowenstein & 
Mushotzky 1996). Similarly, Baldry & Glazebrook (2003) 
find that a slightly flatter IMF slope of F = 2. 15 is favored to 
fit the star formation history (SFH) and luminosity density in 
low redshift star-forming galaxies. Other evidence for a non- 
Salpeter IMF is shown by varying the IMF slope or extending 
the IMF upper-mass limit up to at least 120 M© over the typ- 
ical values of 100 M© in order to reproduce the observed Ha 
equivalent widths in Sloan Digital Sky Survey (SDSS) galax- 
ies (Hoversten & Glazebrook 2008; 2010). Additionally, large 
Lya equivalent widths in some Lya emitting galaxies have 
been observed which cannot be explained by normal stellar 
populations, but may be caused by a top heavy (i.e., flatter) 
IMF slope (e.g., Kudritzki et al. 2000; Malhotra & Rhoads 
2002; Finkelstein et al. 2007). 

One way to study the IMF is to compare SFR indicators 
based on nebular lines, whose luminosities are proportional 
to the number of ionizing photons generated primarily by O- 
stars, to those based on L/«. In a dusty starburst most of the 
ionizing radiation is absorbed by dust and reprocessed into the 
thermal infrared, characterized by L/r, where L/r is expected 
to be proportional to the bolometric luminosity (Lhof, Kenni- 
cutt 1998). While Herschel can permit direct measurements 
of L/R, the majority of galaxies at cosmological distances fall 
well below the 5cr confusion limit of 25-35 mJy for Herschel 
at 250-500 /im (e.g. Chapman et al. 2002; Kneib et al. 2004). 
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TABLE 1 

SPIRE Flux Density Measurements for Component B° 



Wavelength Flux Density Measured EiTor Confusion Error'' Total Quadrature-Summed Error"^ 

[fj.m\ [mJy] [mJy] [mjy] [mjy] 



250 55.5 2 5.8 6.1 

350 56 7 6.5 9.4 

500 34 8 7.3 10.5 



Note. — "Photometry for Component A is approximately 1.5 times fainter. 'Nguyen et al. (2010). '^These total errors do not include an additional 15% error 
based on SPIRE's estimated absolute calibration uncertainty (Griffin et al. 2010). 



Fortunately, gravitationally lensed systems can easily be ob- 
served in the sub-mm by Herschel and other sub-mm tele- 
scopes. One such galaxy is SMM J163554.2H-661225 at z = 
2.515, which is gravitationally lensed by the rich cluster Abell 
2218 (Kneib et al. 2004). The system was first discovered in 
the sub-mm at 450 and 850 /im by SCUBA on the James Clerk 
Maxwell Telescope (JCMT; Kneib et al. 2004) as a multiply 
lensed system. The brightest component (component B) has a 
lensing magnification of 22 and the second brightest compo- 
nent (component A) has a lensing magnification of 14 (Kneib 
et al. 2004). SMM J 1 63554.2H-66 1 225 was observed by Rigby 
et al. (2008) in the mid-IR with Spitzer from 3.6-8.0 /^m with 
the Infrared AiTay Camera (IRAC; Fazio et al. 2004), at 24 - 
160 /im with the Multiband Imaging Photometry for Spitzer, 
(MIPS; Rieke et al. 2004), and with the Spitzer Infrared Spec- 
trograph (IRS; Houck et al 2004). Papovich et al. (2009, here- 
after P09) also analyzed mid-IR Spitzer IRS spectroscopy of 
component B to make the first detection of Paa emission in a 
high-redshift galaxy, from which they derived a revised dust- 
corrected Paa-SFR of 225 ± 37 M© yr"' for a Salpeter IMF 
with an upper and lower mass cutoff of 100 Mq and 0. 1 Mq. 

In this paper, we present new Herschel Spectral and Pho- 
tometric Imaging Receiver (SPIRE; Griffin et al. 2010) ob- 
servations of SMM J163554.2H-661225. In §2 we analyze 
the SPIRE data, including the measurement of crowding- 
corrected flux densities. In §3 we derive the rest-frame FIR 
luminosity, dust temperature, and IR SFR. We compare these 
results to local samples and other high-redshift sub-mm galax- 
ies. In §4 we discuss the implications of the variations in SFRs 
on the star formation history and IMF. In §5 we present our 
summary and conclusions. Where applicable, we use a cos- 
mology with H() = 70 km s"' Mpc"\ 0.^ = 0.3 and fi^ = 0.7. 

2. DATA ANALYSIS 

Herschel SPIRE 250, 350, & 500 pm data of Abell 2218 
were taken as part of the Herschel Multi-tiered Extragalac- 
tic Survey (HerMES; Oliver et al. 2010). The reduced 
SPIRE data of Abell 2218 were accessed through the HeDaM 
databasfl and were processed using the Herschel Interactive 
Pipeline Environment (HIPE v2.3) as described in Oliver et al. 
(2010). The SPIRE 250, 350, and 500 images are blended 
for the two components (A and B, see Figure 1) of the lensed 
galaxy SMM J163554.2H-661225. We used the galaxy fitting 
software GALFIT (v3.0 Peng et al. 2010) to perform point 
source function (PSF) fitting of the SPIRE data for the two 
components of the galaxy. Figure 1 shows a 215" x 215" re- 
gion around SMM J163554.2H-661225 in Spitzer MIPS 24 
p,m, and the SPIRE 250, 350, and 500 pm data, with 250 
/im contours overlaid on the other images. We used the po- 
sitions of the two components from the MIPS 24 /im image 
as a prior for the location of the components in the SPIRE 
250 /im image. We ran Source Extractor (Bertin & Arnouts 



1996) on the SPIRE images prior to GALFIT to estimate in- 
put magnitudes, radial profile, and positions of other sources 
in the image. GALFIT was run on each of the SPIRE im- 
ages on a 2 15" X 2 15" region centered on the two dominant 
lensed components, fitting all sources. GALFIT requires both 
an uncertainty image and a PSF. The uncertainty images from 
the Herschel archive for this field were used, and we adopted 
the HerschellS,?mE theoretical PSF^ The SPIRE 250, 350, 
and 500 /im beams are Gaussian with full width half max- 
imum (FWHM) values of 18.1", 25.2", and 36.6", respec- 
tively. We also compared the results using model PSFs from 
the GOODSI Herschel survey (PI D. Elbaz). The flux density 
measurements varied by less than 5% based on the choice of 
PSF 

The components of the lensed galaxy are reasonably re- 
solved in the 250 /im image (see figure 1). Therefore, we 
measured the flux density for each source in the 250 /im im- 
age using GALFIT, and allowing the position to vary within 
the 1(7 positional uncertainty from SExtractor with no other 
constraints. We measure a flux density at 250 /im for Com- 
ponent B (the brightest component), S250 = 55.5 ± 2 mJy. 
The measured flux density ratio between the two compo- 
nents is S25o(B)/S25o(A) = 1.5 ± 0.1. This ratio is consis- 
tent with the flux density ratio measured at 24 /im from P09, 
S24(B)/S24(A)=1.6, where the 24 /im flux densities of compo- 
nents A and B are 0.72 mJy and 1.16 mJy, respectively. This 
flux density ratio is also consistent with the reported ratios of 
components B and A at 450 and 850 pm (Kneib et al. 2004). 

The A and B components of this lensed galaxy are more 
severely blended in the longer wavelength SPIRE 350 and 500 
/im data (see figure 1). Therefore, we add a priori constraints 
on both the astrometric positions and flux-density ratio to ex- 
tract the individual flux densities. We measured flux densities 
for the components by fixing the positions of each component 
to have the value measured in the 250 pm image. We also 
include a constraint that the flux density ratio at 350 and 500 
/im for the two components have the same ratio measured at 
250 pm. Given the uniformity in the flux-density ratios of 
the A and B components from the mid-IR (24 /im), far-IR 
(250 /im) and sub-mm (450-850 pm) (see paragraph above), 
this assumption seems valid. Including these constraints with 
GALFIT, we derived flux densities at 350 and 500 /im for 
component B of Ssso = 56 ± 7 mJy, and S500 = 34 ± 8 mJy. 
The values for component A are a factor of 1 .5 fainter (by con- 
struction). The SPIRE flux density measurements for compo- 
nent B are detailed in Table 1 . 

In addition to the measured flux density uncertainties, 
SPIRE data suffer noise from confusion, for which Nguyen et 
al. (2010) report Icr errors of 5.8, 6.5, and 7.3 mJy per beam 

' http://hedam.oamp.fr 

^ http://herschel.esac.esa.int/twild/bin/view/Public/SpireCalibrationWeb 
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Fig. 1. — Mid and Far-IR images of the lensed galaxy SMM 
J163554.2+661225. The four panels show 215" X 215" regions around the 
two dominant lensed components (A and B, as labeled) in the Spitzer MIPS 
24, Herschel SPIRE 250, 350, and 500 ^lm data (as labeled). The contours 
show levels 1, 3, and 5 cr above the noise per pixel in the 250 /xm data. While 
the 24/jm data resolve components A and B of SMM 1163554.1+661225, 
these components are blended in the Herschel SPIRE data. Accordingly, we 
fit the Herschel data with GALFIT using the positions of components A and 
B from the 24/jm image. 

for 250, 350, and 500 jivn.. Therefore, we adopt conserva- 
tive errors, adding the flux measurement and confusion errors 
in quadrature, giving total flux density errors of 6.1, 9.4, and 
10.5 mJy. We use these total errors when fitting the suite of 
models and computing IR luminosities and dust temperatures 
as detailed below. We also use previously published MIPS 
70, SCUBA 450 and 850 //m flux densities, which for Com- 
ponent B are 2.56 ± 0.9 mJy, 75 ± 15 mJy, and 17 ± 2 mJy, 
respectively (P09; Kneib et al. 2004). 

3. RESULTS 

3.1. Infrared Luminosity and Dust Properties 

We use the Herschel SPIRE flux densities, Spitzer MIPS 70 
/im data, and SCUBA 450 and 850 /im data to construct an 
infrai-ed SED for Component B of SMM J163554.2+661225. 
We fit a suite of IR SED templates from Chary & Elbaz (2001; 
hereafter CEOl) and Dale & Helou (2002; hereafter DH02) to 
the combined data, finding the best fit model via minimiza- 
tion. The total L/r is derived by integrating the best-fit models 
between 8 to 1000 /im. From the best-fit DH02 template we 
derive a total infrared luminosity of 6.8 x 10" L0, with a re- 
duced = 1-4. For the best-fit CEOl template we derive a 
total IR luminosity of L/r = 7.0 x lO" Lq, with a reduced 
= 1.2. Both values of L/s are corrected for the gravitational- 
lensing magnification factor of 22. Figure 2 shows the IR SED 
and the best-fit SED curves to the data. Figure 2 also shows 
the best-fit to the 70, 450, and 850 /im data alone. Without the 
Herschel data the best fit template has a peak that is shifted 
to longer wavelengths. To determine the uncertainty on our 
best-fit models we ran 10'* Monte Carlo simulations, varying 
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Fig. 1.— The infrared SED of SMM J163554.2+661225. The top panel 
shows flux densities from MIPS 24 and 70 ^tm (triangles), SPIRE 250, 350, 
and 500 /im (diamonds), and SCUBA 450 and 850 ^m (squares) measured 
flux densities. All error bars on the data points are shown at Icr. The solid and 
dashed curves show the IR SED template fits to the 70, 250, 350, 450, 500, 
and 850 /jm flux densities, using templates from DH02 (solid line) and CEOl 
(dashed line). The Herschel data constrains the peak of the dust emission to 
shorter wavelengths compared to the best-fit CEOl template to the SCUBA 
and MIPS 70 [im data alone, shown by the dotted line. The bottom panel 
shows the ratio of the best-fit models and data points to the best-fit model of 
the CEOl model (to all data points). As with the top panel, the error bars are 
Icr. 

the input flux density measurements by a Gaussian random 
amount proportional to their errors in each simulation. The 
resulting 68% confidence ranges of the derived total IR lu- 
minosity from the best-fits are: 6.8 ± 0.6 xlO" and 7.0 
± 0.5 xlO" L0. Combining the results from the different 
templates we get an average L/r = 6.9 ± 0.6 x 10" L©; thus 
the uncertainty on L/r from the photometric data uncertainties 
dominate over systematic effects from the choice of models. 
Previous estimates of L/r ranged from 5 -10x10" L0 based 
on DH02, CEOl, and Rieke et al. 2009 templates (P09) and 
5.7 - 9.5x10" Lq based on DH02 templates (Rigby et al. 
2008). Thus the addition of the Herschel data results in higher 
accuracy, and shrinks the uncertainty on L/r from a factor of 
order 2 down to 10%. 

The CEOl SED templates were constructed to match cor- 
relations between the observed mid-lR and FIR fluxes of lo- 
cal galaxies. The shape of the FIR SED of the CEOl mod- 
els are calibrated on the local IR-luminosity - dust tem- 
perature relationship. As such, the CEOl templates with 
higher dust temperature correspond to galaxies with higher 
L/R. However, we find that the shape of the FIR SED of 
SMM J163554. 2-1-661225 corresponds to a much cooler dust 
temperature compared to local galaxies of similar luminos- 
ity. Indeed, the CEOl template that best matches the FIR SED 
of SMM J163554.2H-661225 has an intrinsic luminosity that 
must be "scaled-up" by a factor of 7.5. Therefore, this best-fit 
template must be scaled-up by a factor of nearly an order of 
magnitude to match both the shape of the SED and the flux 
density measurements of SMM J163554.2H-661225. Similar 
results were seen by Muzzin et al. (2010) for their sample of 
two z ^ 2 galaxies, where their galaxies were comparable to 
local luminous FIR galaxies, but scaled up in luminosity by 
more than an order of magnitude. These new results confirm 
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Fig . 3 . — Total IR luminosity versus dust temperature at z = and z ~ 2. 
SMM J163554. 2+661225 is shown by the black square. The error bars shown 
for SMM J163554.2+661225 are Icr. Black diamonds are z > 1 Herschel- 
ATLAS sources (Amblai'd et al. 2010); blue X's: z > 1 BLAST sources 
(Dye et al. 2009); green asterisks: z ~ 2 SMGs (Chapman et al. 2005); cyan 
crosses: z ~ 2 ULIRGs with Herschel detections (Magdis et al. 2010); red tri- 
angles: local IRAS-selected galaxies with SCUBA 850/^m detections (Dunne 
et al. 2000). 

the conclusions of Rigby et al. (2008), and P09 who used mid- 
and far-IR photometry and mid-IR spectroscopy to demon- 
sti-ate that SMM J163554.2+661225 is inconsistent with the 
spectra of local ultra luminous infrared galaxies (ULIRGs), 
but is consistent with the SED shape of 1ow-L/r local starburst 
galaxies and has just been "scaled-up" in luminosity and star 
formation by 1-2 orders of magnitude. 

Using the best-fit SEDs we calculated an effective dust tem- 
perature. To = 36 ± 3 (Ict) K, assuming a modified grey -body 
distribution (e.g. Young et al. 1989; Calzetti et al. 2000; Am- 
blard et al. 2010) with a fixed emissivity parameter, (3=1.5, of 
the form: 

The uncertainty on the dust temperature was derived by fit- 
ting all of the DH02 and CEOl models within Icr of the best- 
fit, based on the Monte Carlo simulations described above, 
to Equation 1. Based on the mid-IR and SCUBA sub-mm 
data, P09 estimated To = 52 K. We favor the dust temper- 
ature measured here using the Herschel data, as these data 
now constrain the peak of the IR emission, providing the most 
constraining power on the dust temperature. The best-fit dust 
temperature again demonstrates that the SED shape of SMM 
J163554. 2-1-661225 is not consistent with local galaxies of 
comparable L/r, which have significant contributions of warm 
(> 70 K) dust to the IR emission (e.g., CEOl; DH02; Rieke et 
al. 2009). 

We compare the calculated dust temperature and IR lumi- 
nosity of SMM J163554.2H-661225 to other galaxies that have 
measured dust temperatures using Herschel data or other sub- 
mm observations. Figure 3 shows To vs. L/r for galaxies at 
similar redshifts (1 < z < 3), including z > 1 sources selected 
from the He rschel-ALTLAS survey (Amblard et al. 2010), z 
> 1 sources detected in BLAST (Dye et al. 2009), z - 2 
sub-millimeter galaxies (SMGs; Chapman et al. 2005), and 
Herschel detected z 2 ULIRGs (Magdis et al. 2010). We 
also compare to a sample of local IRAS-selected galaxies with 



SCUBA 850^m detections (Dunne et al. 2000). The other 
samples compute L/r (from 8-1000 /im) and To in similar 
ways with a fixed (3 = 1.5, except for the local IRAS sources 
which have varying values of /3. The IRAS sources also quote 
Lf/R instead of L/r, and we converted these using L/r =1.4 
X L/r/R (Dale et al. 2001). If we change (3 over a range of 2.0 
to 1 .2 for our object, then this corresponds to a range in dust 
temperatures of 25 - 42 K. However, if /3 is allowed to vary, 
the best-fit value of /3 = 1.6 results in To = 34.5 K; consistent 
with our result for a fixed /3 = 1.5. 

In Figure 3, SMM J163554.2H-661225 has a similar dust 
temperature to the local IRAS sources, almost all of which 
have lower L/r. It has a dust temperature that is lower 
than local IRAS sources of comparable L/r. While SMM 
J 163554.2-1-66 1225 has a somewhat warmer dust temperature 
than other z ^ 1-2 SMGs of comparable L/r; Chapman et 
al. (2005) discuss how the dust temperature of high-redshift 
sub-mm-selected samples are biased to cooler temperatures. 
Furthermore, SMM J163554.2H-661225 appears to to be in a 
"bridge" region between IR-luminous galaxies with cool and 
warm dust temperatures. Magdis et al. (2010) also note that 
their z ^ 2 ULIRG sample with Herschel detections bridges 
the 'cooler' high-z SMGs to the 'warmer' local/intermediate- 
z ULIRGs. 

We estimate the total dust mass using our estimate of the 
dust temperature and the FIR luminosity. To estimate the total 
dust mass in SMM J163554.2H-661225, we adopt the formu- 
lation of Young et al. (1989) and Calzetti et al. (2000): 

M^,,„ = C5ioo£>'[exp(143.88/W)-l] [Mq] (2) 

Where S loo is the rest-frame flux at 100 fim in Jy, D is the dis- 
tance to the galaxy in Mpc, the expression in square brackets 
is the temperature-dependent part of the blackbody distribu- 
tion, and C is a scale factor that depends on the physical prop- 
erties of the dust grains: C ^ 6 M© Jy"' Mpc"^ for (3 = 2 and 
-5 Mq Jy-i Mpc"- for (3=1.5 (Young et al. 1989; Calzetti 
et al. 2000). For the range in single component dust tempera- 
tures, we estimate a total dust mass for this galaxy of 2.8 g 
(Icr) X 10^ M©, corrected for the gravitational-lensing mag- 
nification factor of 22. This is approximately a factor of 10 
higher than the estimated dust mass from Kneib et al. (2005) 
based on the 850 /im flux density and a dust temperature of 
~ 50 K. This can be explained by the colder dust tempera- 
ture determined here; as there will be a larger fraction of cold 
dust for this galaxy with a dust temperature of 36 K compared 
to 50 K, resulting in a higher total dust mass. Based on the 
molecular gas estimate of 4.5 x 10^ Mq by Kneib et al. (2005) 
this implies a dust-to-gas mass ratio of 0.06. This is consis- 
tent with ratios found for local star forming IR galaxies in the 
Survey of Nearby Infrared Galaxies (SINGS; Kennicutt et al. 
2003; Draine et al. 2007) sample. 

P09 also noted that SMM J163554.2H-661225 has a 
lower ratio of L(24 /im)/L(PaQ!) compared to what is mea- 
sured for any local ULIRGs, demonstrating that SMM 
J163554.2H-661225 lacks a warm (To 70 K) dust com- 
ponent. However, there could be a colder dust component 
present in SMM J163554.2H-661225. The local samples of 
IRAS selected galaxies (Dunne et al. 2000) preferentially de- 
tected galaxies with large amounts of To > 40 K, but because 
these galaxies were selected with IRAS they miss most of the 
cold dust dominated sources. The observations of optically 
selected galaxies of the SCUBA Local Universe Galaxy Sur- 
vey (SLUGS; Vlahakis et al. 2005), supplemented by Spitzer 
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data of a subsample of SLUGS galaxies by Willmer et al. 
(2009), showed that the mid-lo-far-IR SEDs of these galaxies 
are better matched with both cold and warm dust components. 
These samples show on average a 2 x larger proportion of cold 
dust relative to warm dust than solely IR-selected galaxy sam- 
ples do, such as the IRAS selected sample (Dunne et al. 2000) 
and SINGS sample (Kennicutt et al. 2003). 

We fit models of Draine & Li (2007) to the SED of SMM 
J163554. 2+661225 in order to investigate a two component 
dust fit. The impUed total IR luminosity from the best-fit 
Draine & Li model is 7.0x lO" Lq, consistent with our other 
SED fits. The best-fit Draine & Li (2007) SED is modeled 
with a modified grey-body distribution with two dust compo- 
nents, each with /3=1.5. This fit corresponds to warm and cold 
effective dust temperatures of 46 K and 28 K. Using the dust 
mass equation above, we calculate warm and cold dust mass 
components of Mcu = 4.3 x 10^ Mq, and M„a„„ = 5.6 x 10^ 
Mq. The mass from the warm dust component more closely 
matches the estimated total dust mass determined by Kneib et 
al. (2005). The ratio of Mcoid/^warm = 7.8, with a total dust 
mass = 4.9x10^ M0, consistent within 2a of our result of 
the dust mass estimate from the single-component tempera- 
ture fits. This lower ratio of cold dust to warm dust for SMM 
J163554.2+661225 is more similar to that of the SINGS sam- 
ple than in the SLUGS sample (Draine et al. 2007; Willmer 
et al. 2009). The SLUGS sample typically had higher ra- 
tios of McoidM„arm, ou Order -1000 (Willmer et al. 2009). 
This agrees with the result of Rigby et al. (2008) who noted 
the similarity between the mid-IR emission features of SMM 
J163554.2+661225 and that of local IR star-forming galaxies. 

4. STAR FORMATION RATES AND IMPLICATIONS FOR IMF 

Local star-forming galaxies show a tight trend between 
L(Paa) and L/« (e.g., Calzetti et al. 2005; Alonso-Herrero 
et al. 2006), both of which trace the instantaneous SFR. We 
calculate the SFR of SMM J163554.2+661225 based on the 
derived total IR luminosities from the SED fitting, using the 
relation derived by Kennicutt (1998), which gives SFR = 
4.5 X lO""*"* X FIR (erg s"') = 119 ± 10 Mq yr"' for a Salpeter 
IMF from 0.1 - 100 M©. Previous estimates of SFRs from 
the total IR emission in this galaxy range from SFR/r = 90 - 
180 Mq yr-1 (P09) and SFR/r = 140 ± 30 M© yr"' (Rigby 
et al. 2008). Our estimate of the IR emission from Spitzer, 
Herschel, and the sub-mm are consistent with these results, 
although the uncertainty on our measurement is significantly 
reduced by including the Herschel far-IR data. 

P09 obtained a Paa measurement of component B using 
Spitzer IRS spectroscopy with the SL2 module and measured 
a Paa-SFR; the Spitzer IRS observations only covered com- 
ponent B, so we can only compare the IR-derived SFR to the 
Paa-SFR for component B. P09 measured a Paa line flux for 
component B of 8.6 ± 1.4 x 10"'^ erg s"' cm"^, which is un- 
corrected for dust attenuation and gravitational-lensing mag- 
nification. Kneib et al. (2004) measured an Ha line flux of 
5.9 X 10"''' erg s"' cm"^. We re-calculate the dust-corrected 
Paa luminosity, using an updated value for the dust attenua- 
tion for the Paa Une (see P09). We use the ratio of the Paa 
line flux to the Ha line flux to estimate the amount of dust 
attenuation affecting the nebular gas, assuming the Calzetti et 
al. (2000) dust law. This results in an extinction estimate of 
A(V) = 3.9 ± 0.4 mag, and this corresponds to A(Paa) = 0.58 
± 0.04 mag (we assume that the stellar continuum is attenu- 
ated at the same level as the nebular lines [e.g. Erb et al. 2006; 



though see also Forster Schreiber et al. 2009]). Correcting for 
the dust extinction, and applying a gravitational-lensing mag- 
nification factor of 22, the Paa Une luminosity is L(Paa)cor = 
3.38 ± 0.55 X 10''^ erg s~^ This corresponds to an extinction 
and lensing-corrected SFRpa„ = 225 ± 37 Mq yr"', following 
the Kennicutt (1998) relations. This luminosity corresponds 
to an ionizing continuum flux of Q,o„ = 2.1 ± 0.3 x 10^^7 
s~^ This updated dust correction results in a 30% increase in 
the derived Paa-luminosity and Paa-SFR, and Q,„„ from that 
originally derived by P09. 

Our estimate of the SFR derived from the IR emission, 
SFR(IR) = 1 19 ± 10 Mq yr"' is significantly different than the 
SFR calculated from the dust-corrected Paa line, SFR(Paa) = 
225 ± 37 Mq yr"i. They are offset by a factor of 1 .9 ± 0.35, a 
2.5(7 discrepancy (or a 99% confidence). The offset between 
the SFR derived from the Paa line and IR implies that one 
or more of our assumptions are incorrect. We consider here 
several explanations to account for the apparent discrepancy 
between the IR and Paa SFRs. 

4.1. Variations in Dust Attenuation 

In our analysis, we corrected the Paa line luminosity as- 
suming a Calzetti dust extinction law (Calzetti et al. 2000). 
The exact level of dust attenuation is uncertain (see below), 
but changes in the assumed extinction are unable to account 
fully for the measured offset between the IR and Paa SFRs. 
For example, P09 report a Paa luminosity of (2.05 ±0.33) x 
10"*^ erg s"' uncorrected for dust extinction. This corresponds 
to a SFR of 130 ±21 M© yr"^ If fliere is no dust extinc- 
tion correction the SFRs are consistent at the Icr level. How- 
ever, zero dust attenuation seems unlikely. First, the ratio of 
the measured nebular emission lines are strongly inconsis- 
tent with the assumption of zero dust (see below). Second, 
P09 showed from modeling the rest-frame UV-to-NIR pho- 
tometry of the galaxy with stellar population synthesis mod- 
els (Bruzual & Chariot 2003) that the dominant star-forming 
component has a best fit with A(V) — 3.2 mag. Finally, the 
simple fact that this galaxy is bright in the far-infrared im- 
plies that it contains a significant amount of dust (see §3.1). 
Therefore, significant dust extinction seems required. 

In our analysis, using the dust extinction law of Calzetti et 
al. (2000) results in a line extinction measurement of A(Paa) 
= 0.58 mag derived from the observed ratio of the Ha and 
Paa lines. Using other dust attenuation laws (e.g., Dopita 
et al. 2005; CardelU et al. 1989) results in slightly larger 
extinction estimates, and would therefore increase the dust- 
corrected Paa luminosity by —10%. This amplifies the offset 
between the Paa-derived SFR and IR-derived SFRs. Thus, the 
use of the Calzetti dust law results in a conservative estimate 
of the Paa luminosity and SFRpaa, and therefore minimizes 
the differences between the two SFR measurements. 

The Paa dust extinction determined above is based on the 
ratio of the Ha line flux from Kneib et al. (2004) to the Paa 
line flux. However, there are some discrepancies in the re- 
ported Ha line flux measurement, and it could be as much as 
4x fainter (P09, W. Rujopakarn et al., in preparation), but a 
fainter Ha line flux would image a larger derived A(V) and 
A(Paa), and would increase the extinction-corrected Paa lu- 
minosity and SFR (e.g., a Ha line flux lower by a factor of 4 
would increase dust-corrected Paa line flux by — 30%). 

We also we re-calculate A(Paa) using a newer measure- 
ment of H/3 for this galaxy from Richard et al. (201 1). The ob- 
served Paa line flux for this galaxy again is 8.6 ± 1.4 x 10"'* 
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erg s" , uncorrected for the gravitational lensing magnifica- 
tion or dust attenuation (P09). The observed Ha line flux for 
component B (making no correction for uncertainties in grav- 
itational lensing, dust attenuation, or slit losses) is 5.9 x 10""' 
erg s-i (Kneib et al. 2004), and the H/3 line flux is 6.88 ± 
0.29 xlO"'^ erg s"' (Richard et al. 2011). Given these three 
lines we calculate an average E(B-V) of 0.97 ± 0.07 and this 
corresponds to an average A(Paa) of 0.57 ± 0.06 assuming 
the Calzetti et al. (2000) attenuation law. Therefore, the ad- 
dition of the H/3 line confirms the estimate of the dust atten- 
uation affecting the nebular gas, and implies the extinction is 
approximately optically thin to these photons. 

Lastly, the extinction correction derived here assumes the 
nebular gas is approximately optically thin to photons from 
Ha and Paa. While this is supported by the measured line 
ratios, if the opacity of the nebular gas is substantially higher, 
then it will obscure a higher fraction of emission associated 
with star formation. However, if this is the case, then the 
A(Paa) we derive above is only lower limit, as it is only sen- 
sitive to the outermost layer of the gas. If this is the case, then 
the dust corrected Paa luminosity will also be higher, which 
further exacerbates the problem. 

Therefore, in summary it seems unlikely that assumptions 
about the dust law or the measurements of the dust attenuation 
of the Paa line contribute significantly to the observed offset 
between the IR and Paa SFRs. Indeed, our assumptions about 
the dust attenuation are mostly conservative and the intrinsic 
Paa may be larger than reported here. 

4.2. Variations to the ISM Conditions 

The calculations above depend on intrinsic line ratios as- 
suming Case B recombination at T^, = 10** K. We also con- 
sider how changes in the assumed ISM temperatures effect 
the estimation of A(Paa), and therefore the derived Paa lu- 
minosity. Increasing the ISM temperature to larger values (T^, 
= 2x 10"* K), as could be expected for a highly star-forming 
galaxy, increases the intrinsic value of Ha / Paa from 8.46 to 
9.68 (Osterbrock 1989). This increases the estimated value of 
A(Paa) to 0.62 mag, an increase of 10%. This would slightly 
increase the dust-corrected values of L(Paa) and SFRfaa- 

As an extreme case, A(Paa) would be reduced to 0.46 mag 
if one assumes Case A recombination and a very low ISM 
temperature of T^ = 2500 K. This would decrease L(Paa)co,- 
and SFRpaa by a maximum of 12%, but it is insufficient to 
account for the factor of ^2 difference between the Paa and 
IR SFRs. Therefore, our assumption of Case B recombina- 
tion with an ISM temperature of 10** K does not affect our 
conclusions. 

4.3. Changes to the Star Formation History 

The offset between the two SFRs could also be due 
to assumptions about the stellar population of SMM 
J163554.2H-661225. The Paa and SFR calibrations both 
assume a constant SFR and a stellar population age of 100 
Myr, as well as a Salpeter IMF. Using the 2007 version of the 
Bruzual & Chariot (2003) stellar population synthesis mod- 
els, we investigated if this factor of 1.9 between the derived 
SFRs could be accounted for by varying the star formation 
history or age, but still with a Salpeter IMF. For example, 
the commonly assumed relations of Kennicutt 1998 assume a 
roughly constant star-formation history for the past 100 Myr 
However, there is evidence that galaxies at z > 2 have star- 
formation histories that increase with time (e.g., Papovich et 



al. 2011), which affects the relative proportion of stars that 
contribute to nebular emission and the IR emission. 

From the stellar population synthesis models with a 
Salpeter IMF, we calculate the ratio of ionizing photons (Q,o„) 
to bolometric luminosity versus stellar population age. We as- 
sume that the ratio of Q,,,,, / Li,,/ is proportional to the ratio of 
Paa luminosity to total IR luminosity (as argued by Kenni- 
cutt 1998). We tested how a SFR that rises with time affects 
the ratio of Q,o„ / 'Lboi- Assuming the SFR rises with time 
(approximated by the delayed SFR model within Bruzual & 
Chariot 2003), we find a maximal increase of Q,,,,, / hi,oi by at 
most a factor of 1.25 compared to this ratio for constant SFR. 
Other increasing SFRs (including a model with exponentially 
increasing SFR) were also investigated which produce more 
Ly man-continuum photons by as much as a factor of 1.3 for 
ages < a few e-folding times. Higher ratios are possible, but 
only in extreme situations, such as exponentially increasing 
SFRs at ages >> a few e-folding times, but it seems physi- 
cally unlikely that a galaxy will increase its SFR exponentially 
for such sustained periods without disruption (i.e. feedback). 
Therefore, star formation history alone can only partially ex- 
plain the offset between the two derived SFRs. 

4.4. Varying the IMF 

One possible way to create more ionizing photons would 
be an IMF that is weighted towards high-mass stars compared 
to that of an IMF with a Salpeter-like IMF slope. We inves- 
tigated how changes in the the form of the IMF affects the 
Q,o„ / hhoi ratio, including an IMF with a flatter slope (i.e., a 
"top-heavy" IMF) or an IMF with a higher upper mass cut- 
off. The expected ratio of the total number of ionizing pho- 
tons (Q,o«) to the total IR luminosity, Q,o„ / L/^, following 
the Kennicutt (1998) relations and assuming a Salpeter IMF 
from 0.1 to 100 Mq, should have a ratio of 1.6 xlO'*-' 7 s"' 
L0"'. This assumes that the Paa-derived SFR and the IR- 
derived SFR should be equal. However, our measured ratio 
for SMM J163554.2H-661225 is Q,„„ / Ur is 3.0 ± 0.6 x lO'*^ 
7 s"' L0"', reflecting the factor of 1.9 difference between the 
two SFRs compared to the theoretical value. 

We computed Q,,,,, and L/r (assumed to be ^ Li,,;) as a 
function of both the IMF slope and upper-mass cutoff. To 
make this calculation, we take empirical measures of Qo as a 
function of spectral type for 0-stars with masses 19.3 - 100 
M0 derived by Sternberg et al. (200313. We then weighted 
the values by the number of stars per unit mass for each IMF, 
and integrate to calculate the total number of ionizing photons 
for a given IMF. We repeated the calculation for the bolomet- 
ric luminosity, assuming a mass-luminosity relation of L cx 
M^ for stars less than 20 Mq (Popper 1980), and L cx M' 
for M > 20 Mq based on empirical models of OB stars (Stern- 
berg et al. 2003). 

Figure 4 shows the effect on the ratio Q/o,, / L/« for IMFs 
of different forms, including changes in the IMF slope and 
the upper-mass cutoff. We have overplotted both the expected 
ratio from Kennicutt (1998) and our measured ratio. The left 
panel shows that a shallower IMF slope, with a best-fit of F = 
1.9 ± 0.15, reproduces our measured ratio of Q,o„ / L/r. This 
result is a 3cr variation from Salpeter (F = 2.35). As shown 
in the right panel, a Salpeter IMF slope with an upper-mass 
cutoff extended up to ^200 Mq also reproduces the measured 

' We assume that Qo increases with stellar mass at a rate derived from the 
empiiical measurements of Sternberg et al. (2003). We extrapolated this rate 
to estimate Qq values above 100 Mq. 
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Fig. 4. — Effect of changes to the IMF on the ratio of the number of ionizing photons (Qum) to L//;. Both scenarios or a combination can reproduce the 
measured difference between the SFR;r and SFR/j„q ■ The left plot shows how varying the IMF slope changes the derived ratio of Q,o„ to L//;. The red dot-dash 
line shows the expected ratio following the Kennicutt relations assuming a Salpeter mass distribution with upper and lower mass cutoffs of 100 Mq and 0.1 Mq. 
The blue dashed line shows our measured value of Q,,,,, / L/r, the blue shaded area shows the Icr uncertainty on the measured ratio. To match our observations, 
more ionizing photons need to be made relative to the L/r (traced by the bolometric luminosity, see Kennicutt 1998) which can happen with a shallower IMF 
distribution with slope F = 1.9 ± 0.15. The right plot shows that we can also reproduce our observed ratio of Q,o„ / L/r with a Salpeter IMF by extending the 
upper-mass limit to 200 Mq. The grey hatched region in each panel shows our preferred Icr values. 



ratio. The la uncertainty on the measured ratio of Q,o„ / L/s 
is also shown in Figure 4 by the blue shaded area. The vertical 
hatched region in each panel shows our preferred la values. 
Therefore, a flattening of the IMF slope to F = 1.9 ± 0.15 or 
an increase in the IMF upper-mass cutoff to at least 160 M© 
(or a combination of both) is able to account for the observed 
Paa line luminosity and IR luminosity. 

P09 estimate that SMM J163554.2H-661225 has a total stel- 
lar mass of ^ 10'" Mq. This implies that the total number of 
O stars (M > 20 M©) is about 10^, assuming a Salpeter IMF 
from 0.1 to 100 Mq and a constant SFR. If the IMF slope 
varies from the nominal value of Salpeter to F = 1.9 then the 
number of O stars between 20 -100 M© would increase by a 
factor of 3.5. Alternatively, this additional number of O stars 
can be accounted by extending the upper mass cut off to at 
least 160 M©, this would result in ^5% more O stars above 
100 Mq, in which case the majority of the additional ioniza- 
tion comes from these highest mass stars. 

4.5. Metallicity Effects 

We also investigate if changes in metallicity can account 
for the measured discrepancy in the SFRs. In the above cal- 
culations for a varying SFH or a changing IMF, we assume 
solar metallicity. Using stellar population synthesis models 
from Starburst 99 (Leitherer et al. 1999), we specifically try 
to determine if lower metallicities can reproduce the observed 
ratio of Q,o„ / Lto/. We ran four simulations in Starburst 99, 
only varying the metallicity. In all four simulations, we as- 
sumed a constant SFR of 120 M© yr~', and a Salpeter IMF 
between 0. 1 - 100 M©. The four simulations had metallicities 
of 0.02 Zq, 0.2 Zq, 0.4 Zq, and solar metallicity. Based on 



these results, we found that the simulation with solar metal- 
licity at an age of 100 Myr reproduces the expected ratio of 
the number ionizing photons to hb„i (1.6x10"*^ ionizing 7 s~' 
L©"') as calculated following the Kennicutt (1998) relations 
as described above. This is the factor of 1.9 lower than our 
measured ratio of Q,o„ / l^toi- In order to reproduce our mea- 
sured ratio of Q,o„ / l^boi, the metallicity must be lowered to < 
0.02 Zq. Therefore, a very low metallicity population could 
reproduce the higher rate of ionizing photons as compared to 
the total luminosity without changing the IMF. 

However, we have reason to believe that the metallicity of 
this galaxy is nearly solar Using the line ratios of [N II]/Ha 
= 0.3 ± 0.1 (Kneib et al. 2004) and [O III] /H^ = 1.5 ± 0.11 
(Richard et al. 201 1) we calculate the metallicity based on the 
N2 and 03N2 indices. From the N2 index, we get a value of 
12 -I- log(0/H) = 8.6 ± 0.1, assuming the relation from Pettini 
& Pagel (2004), consistent with the solar value 12 + log(0/H) 
= 8.66. From the 03N2 index, we get a value of 8.5 ± 0.06, 
again close to solar metallicity. Therefore the gas metallic- 
ity of SMM J163554.2H-661225 appears to be consistent with 
solar metallicity, and is inconsistent with a metallicity of 0.02 
Z© by more than Aa. Since they have just formed, we can also 
expect the massive OB stars that produce the increased frac- 
tion of ionizing photons to have a similar metallicity as the 
gas, and therefore it is not likely that a very low metallicity 
stellar population is responsible for the high measured value 
of Q,o,i / Lfco/. 

4.6. Differential Magnification 

As a final possibility for a cause of the differing SFRs we 
look at the lensing model of SMM J163554.2H-661225 in or- 
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der to determine if there is evidence for differential magni- 
fication between the physical region of the galaxy that pro- 
duces the Paa emission relative to the IR producing regions. 
For example, if Paa and IR emitting regions have different 
gravitational-lensing magnification factors, then it could con- 
tribute to the discrepancies we observe. 

The gravitational lensing model at the astrometric position 
of the galaxy shows variations in magnification of order 10- 
20% over the 1 arcsecond scale of the galaxy (J.-R Kneib, 
private communication, 2011), and this could contribute to 
the measured offset. However, it seems more likely that the 
Paa and IR emitting regions are coaligned. We note that the 
Spitzer/IRAC data and MIPS 24 fim centroids are aligned to 
better than 1 arcsecond, suggesting the regions contributing 
to the emission at these wavelengths are magnified equally. 
Therefore, while we can not exclude that variations in gravi- 
tational lensing contribute, we argue that it is unlikely to dom- 
inate the measured offset. 

5. SUMMARY AND FINAL THOUGHTS 

The Herschel SPIRE data of SMM J163554.2H-661225 al- 
low for a highly accurate fit to the infrared SED of this 
z=2.515 lensed galaxy. We measured a best-fit IR luminosity 
of 6.9 ± 0.6 X 10" L0, corresponding to a dust temperature 
of 36 ± 3 K, and a dust mass of -3 x 10** Mq. The IR-derived 
SFR is 1 19 ± 10 M© yr"' which is a factor of 1.9 lower than 
the previously derived Paa-SFR. In order to account for the 
discrepancy in the derived SFRs, we have analyzed several 
different scenarios that could alter the ratio of Paa luminosity 
and L/R. One possibility is that there is a varying SFH due 
to a rising SFR. Using a SFR that increases with time, with 
a Salpeter IMF can increase the number of ionizing photons 
relative to the total bolometric luminosity (and thereby to the 
L/r), by as much as 25-30% at a given stellar population age, 
but cannot by itself account for the measured difference be- 
tween SFR/s and SFRpoQ. 

Varying the high-mass end of the IMF can explain our re- 
sults, either by increasing the high-mass cutoff up to at least 
160 M©, or changing the IMF slope to F = 1.9 from the 
Salpeter value of F = 2.35, or a combination of the two. Ei- 
ther scenario produces an increase in the number of massive 
stars, resulting in higher rates of ionizing photons. This ap- 
pears necessary to explain the higher Paa-derived SFR over 
the IR-derived SFR. However, with current data we are unable 
to distinguish between these two scenarios. Similar changes 
in the IMF at the high mass end have been shown to have 
an effect on the relative contributions between various SFR 
indicators (e.g., Daddi et al. 2007), whereas extending the 
IMF upper-mass limit up to at least 120 Mq has also been 
shown to be important in reproducing the observed Ha flux 
in star-forming galaxies from the SDSS (Hoversten & Glaze- 
brook 2008; 2010). Based on our analysis, we favor a scenario 
where the IMF may require minor adjustments to the standard 
form of the Salpeter IMF for some star-forming galaxies. 

Flattening of the IMF at higher redshifts has been proposed 
by other studies to account for offsets between stellar mass 
and SFR densities (e.g., Dave 2008 and references therein). 
A flat IMF, with F = 1, has also been proposed to account for 
SMG number counts (Baugh et al. 2005; Lacey et al. 2008). 
Our results of SMM J163554.2H-661225 aflow for an IMF 
slope r = 1.9 ± 0.15, only slightly flatter than Salpeter Addi- 
tionally, it could be that both a slight change in the IMF and a 
change to the SFH are at play. If the SFH is changing or there 
has been a burst in star formation activity, this could mimic 



some of the changes produced via a varying IMF. However, 
a more complex star formation history only cannot explain 
the whole effect we see; some alteration of the IMF is still 
needed. As discussed above, if the galaxy has a SFR that is 
increasing with time, then this would increase the expected 
number of ionization photons by 25% compared to a constant 
SFR. In this case, coupled with the increasing SFR then the 
IMF would only be slightly flatter than Salpeter with a slope 
r ~ 2.06 and/or with an upper-mass cutoff of up to 120 M©. 

We also investigated if changes to the assumed metallicity 
and dust content of this galaxy could account for the offset 
in derived SFRs. This galaxy would need a very low metal- 
licity population of ^ 0.02 Z©, however this low metallic- 
ity does not match with the observed line ratios of the gas, 
which give a measured metallicity of ^ solar A lower dust 
extinction of A(V) =1.0 mag, corresponding to E(B-V) = 
0.25 mag, and A(PaQ;) = 0.15 mag, could also account for the 
discrepancy between the IR-derived and Paa-derived SFRs. 
However changing the assumptions about the ISM conditions 
(ie. Temperature or Case A vs. Case B recombination) only 
have a small effect on the derived A( V) values for this galaxy, 
and at most can only decrease the derived Paa-SFR by 12%, 
again not enough to make up the difference between the two 
SFRs. Also from stellar population synthesis modeling, pre- 
vious authors (P09) have shown that this galaxy is best-fit 
by a two-component stellar population fit, with the dominant 
star-forming component being very dusty, corresponding to 
E(B-V) = 0.8 mag and A(V) ^ 3.2, and therefore inconsis- 
tent with A(V) ^ 1 . We therefore conclude that the assump- 
tions we have made about the metallicity and dust content of 
this galaxy are consistent with the known observations of this 
galaxy, and it would take extreme changes to one or the other 
in order to be able to explain the differences in the derived 
SFRs. 

Another possibility to explain the offset is that the regions 
of the galaxy that emit the Paa and IR may experience differ- 
ent amounts of gravitational lensing magnification. However, 
we expect the Paa emission to trace the star-forming and IR 
dominated regions. Given the size of the galaxy, this effect 
seems to be less than w 20%, and likely is not a dominant 
effect. 

As a final thought, it may be that a combination of all 
the effects (variations in the star-formation history, minor 
adjustments to the IMF, variations in extinction, etc.) may 
all contribute to the difference in the Paa-derived SFR and 
the IR SFR. Future observations of z^2 lensed galaxies that 
have both Paa emission and FIR observations may be able 
to determine if this offset is common at Z'^2 or if SMM 
J163554.2H-661225 is an unusual case. 
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